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We study the diffusion and deformation of classical 
solitons coupled to thermal noise. The diffusion co- 
efficient for kinks in the (f> 4 theory is predicted up to 
the second order in kT. The prediction is verified by 
numerical simulations. Multiskyrmions in the vector 
0(3) sigma model are studied within the same formal- 
ism. Thermal noise results in a diffusion on the mul- 
tisoliton collective coordinate space (moduli space). 
There are entropic forces which tend, for example, to 
bind pairs of solitons into bi-solitonic molecules. 

Solitons or extended objects are an important ingredi- 
ent in many physical phenomena ranging from bioener- 
getics to superconductivity or nonlinear optics, see e.g. 
|jj for a review. Most of these phenomena take place at a 
finite temperature and involve dissipation and noise. Al- 
though it is sometimes possible to study transport prop- 
erties of solitons starting from first principles M , it seems 
that in most cases one has to rely on some kind of an ef- 
fective diffusive nonlinear equation derived by expansion 
in derivatives of the order parameter, see e.g. |J| for an 
example. 

In the following we develop a perturbative expansion 
in powers of temperature, which we use to predict the dif- 
fusion coefficient and the deformation of kinks at finite 
temperatures. Our predictions are confirmed by numer- 
ical simulations. Then we generalize the formalism to 
multisolitons in the planar vector 0(3) sigma model. In 
the absence of thermal noise solitons of the pure sigma 
model do not interact with each other. Thus they pro- 
vide a convenient setting for an "in vitro" study of noise 
induced entropic forces. 



I. DIFFUSION OF KINKS 

Let us consider a dissipative version of the </) 4 theory 
in one spatial dimension, which is defined, in appropri- 
ate dimcnsionless units, by the stochastic nonlinear field 
equation 



Yd t (j) = dl4> + 2[1 - (/) 2 }(t> + v(t, x) 



(1) 



where L is a dissipation coefficient and rj(t, x) is a gaus- 
sian white noise with correlations 
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< T)(t,x) > = , 

< T}(ti,xi)Tl(t2,X2) >= 2kTY5{t 1 -t 2 )5{xi 



*a) ■ (2) 



The system is coupled to an ideal heat bath at temper- 
ature T. At nonzero temperature the field (f>(t, x) per- 
forms a random walk in its configuration space. In the 
absence of noise, at T = 0, Eq.(^j) admits static kink so- 
lutions </>(t, x) — F(x) = tanh(x). Antikinks are given by 
F(-x). 

Spectrum of kink excitations. Small perturba- 
tions around the kink take the form e~ Jt ^ r u(x). Lin- 
earization of Eq.(|l|) with respect to u(x), for i](t, x) = 0, 
gives 



7 u(x) = --7-5 u(x) 



[4- 



G 



cosh (x) 



u(x) 



(3) 



The eigenvalues and eigenstates can be tabulated as 



7> 
0, 

3, 



u(x) 
F'(x) 

B(x) 

u k {x) 



cosh (x) 
_ sinh(.T) 



cosh 2 (a;) 
ik 3zfctanh(x) -3tanh 2 (z) 
[ + 1 + fc 2 J 



(4) 



where A: is a real momentum. The zero mode (7 = 0) is 
separated by a gap from the first excited state (breather 
mode) of 7 = 3. The continuum states are normalized so 
that 



+00 



dx ul.{x)uk> (x) 



2tt S(k - k') = N{k) 8{k - k') 



1 + fc 2 



(5) 



In the following we will often denote the states (jj) by 
u a (x) with the index a = reserved for the zero mode. 

Collective coordinates. The field in the one kink 
sector of the theory (|l|) can be expanded in the complete 
set (^) as 



4>(t, x) = F[x - £(t)] + ]T A a (t) u a [x - £(*)] 

a=£0 



(6) 



Substitution of the above to Eq.([l]) and projection on the 
orthogonal basis @) gives a set of stochastic nonlinear 
differential equations 
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r n q i(t) - r i(t) Mt) M aQ + Mt)A c (t) p bcQ + 

A b (t)A c (t)A d (t) R bcd0 = r)o(t) , 
r N a A a (t) + la N a A a {t) ~ r i{t) A b (t) M ba + 

A b (t)A c (t) P bca + A b (t)A c (t)A d (t) R bcda = Va{t) 

where ' = d/dt and the coefficients are 



dx u a {x)u* a {x) 



which is singular for £—►£'. According to the first of 
Eqs.([ll|), the probability P(t, £), that the kink random 
walks to £ at the time t, satisfies the diffusion equation 



dP _ d 2 P 



(13) 



N a = 
M ab = 

Pabc = 6 
Rabcd = 2 



+oc 



dx u' a (x)ul{x) 



dx F[x) u a (x)u b (x)u*(x) 



dx u a (x)u b (x)u c (x)u* d {x) 



The ?7 a 's result from the projections 

f + OO 



Va(t) = 



dx rj(t, x)u* a (x) 



(8) 



(9) 



With the correlations (Q) and the orthogonality of the 
basis (^) the correlations of the projected noises are 

< r) a (t) >=0 , 

<V* a (h)Vb(t2) >=2kTTN a 5 ab 5(t 1 -t 2 ) . (fO) 

Diffusion coefficient to leading order in kT. 

The natural parameter of expansion at low temperature 
is y/kT. The projected noises, compare Eq.(|fC|), are 
just of this order. Let us make the customary rescal- 
ing y/kT — > ey/kT. Eqs.(Q) can be expanded in pow- 
ers of e and solved in power series of e with e set to 
f at the end of the calculation. With the expansion of 
the collective coordinates, £ = e£W _|_ g 2 £( 2 ) + . . . and 
A a = sA^ + e 2 A a 2 ^ + . . ., the equations (0) become, to 
the leading order in s, 



r n i (1) (t) = i] (t) , 

TN a A^(t)+ 7a N a A£\t)= Va (t) 



(11) 



In this approximation the zero mode £W (t) and the ex- 
cited modes A^ (t) are uncorrelated stochastic processes 
driven, respectively, by their mutually uncorrelated pro- 
jected noises. 

^'(i) is a Markovian Wiener process whose only non- 
vanishing single connected correlation function is 



<e (1) we (1) (0 >= 



2kT 



5(t - t') 



(12) 



with the diffusion coefficient D = kT/TN Q = 3kT/4T. 
■ Cff^'s are Ornstein-Uhlenbeck random noises with re- 
laxation times r/7 . They represent memory effects, ff, 
for t < 0, T = and then the system is put in contact 
with a heat bath of temperature T > 0, for t > the 
correlations of A^'s would grow according to 



<[A£\t)]* A^(t')>=5 ab 



kT , 



-1a\t-t'\ 



Nala 



[1 



-7a(t+t') . 



(14) 



as can be deduced from a formal solution of the second of 
Eqs.(lf]). All other single connected correlations of A^'b 
vanish. A given mode achieves the state of equilibrium 
with the heat bath after its characteristic relaxation time. 
The longest relaxation time is that of the breather mode. 

The next to leading order correction. Further 
terms in the expansion of collective coordinates can be 
recursively worked out as 



PbcQ 
TNn 



Mo a Jo 



E 

b,c=£0 



Pbc 



dre-^A^(r)A^(r) 



^w = EK (1) w4 2) (t) + (i-2)] 



J2[A^(t)A^(t) + (l<-2)] 

b,c^0 

J2 Ai l \ t )AU(t)4Ht) 

b,c,d^0 



M a0 _ 

No 

PbcQ 



TNn 



)M iWtA Rbcd0 

TN 



(15) 



In the thermodynamic equilibrium the excited modes de- 
velop nonzero expectation values 



<4 2) (°°) >= -J^i\p BBa 



+oc 



dk 



k. — k.a 



j(k)N(k)> ' 
(16) 



where the index U B" stands for the breather mode. In 
particular the expectation value of the breather mode 



is negative, < A% ] [aa) >= -kT^-. The kink gets 
thicker due to the noise. At the same time the expec- 
tation value of <f>[x,t), far from the kink, gets smaller, 
< #»,oo)) >= 1 - &f + 0[(kT) 2 }. 
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The equilibrium correlations of £, up to the next to the 
leading order, are 

<mi(t') >«e 2 <i w (t)i^(t')> + 

e A < £ (2) (*)£ (2) (0 > +e 4 [< £ {1) (t)P\t') > +(1 «- 3)] = 
<C (1) (^ (1) (0 > x 
[1 



3fcT v |M a0 | 2 



7a^a 



— 5>f a0 <A( 2 >(oo) >] + 

a^O 



where . . . denotes terms which are regular as t — > t'. In 
the last equality we have already set e = 1 . The integrals 
over the continuous part of the spectrum and the sum- 
mation over the breather mode yielded < £(t)£(t') >= 
S(t - t') ^f[l + kT 1.8164] + 0[(kT) 3 }.. The diffusion 
coefficient in the equilibrium state, understood as an av- 
erage over times much longer than the relaxation time of 
the breather mode, has been found to be 



D(oo) 



3fcT 
IF 



[1 + kT 1.8164] + 0[kT 3 } 



(18) 



Numerical experiment. The projection of the ex- 
pansion (0) onto F' (x — y), 



I(t,y) 



4-oc 



dx <j)(t, x)F'[x — y] 



(19) 



vanishes for y = £(t). At moderate temperatures (small 
A's) I(t,y) > for y > £(*) and I(t,y) < for 
y < £(i).Thus can be defined as the solution of 
= 0- This definition was introduced in Ref. ||; 
the coordinate defined in this way is a position the kink 
would relax to if the noise were suddenly switched off. 

We performed numerical simulations of Eq.([l]) with 
r = 1 for a range of temperatures. Fig.l compares nu- 
merical results with the perturbative prediction ([L8|), 

In principle, for some periods of time, the solution of 
I[t, £(t)] = may happen to be not unique even for ar- 
bitrarily small T > 0. There is always small but finite 
density of thermally activated kink-antikink pairs. Even 
if this happens for a given realization of the stochastic 
process, one can still find two times t\ < t% such that 
the solution is unique and choose at least one solution 
£(t) valid for any t\ < t < t%\ due to topology the equa- 
tion has at least one solution for any t. The ambiguity 
manifests itself in the first of Eqs (J7]) in the coefficient 
of i, namely T [N Q - Ea#o^(i) M o0 ]. The coefficient 
vanishes and leads to ambiguity, iff I'[t, = or a few 
zeros meet at £(i). At such a point, the trajectory splits 
into several branches or a few branches merge into one. 
The perturbative expansion overcomes this problem, it 
gives an asymptotic approximation to the diffusion coef- 
ficient of the chosen line £ (t) . 



II. DIFFUSION OF MULTISOLITONS. 

A generalization of the ideas developed for a kink to a 
soliton in higher dimensions is straightforward. Less triv- 
ial is the situation with many solitons, which can overlap 
or pass through each other. Let us consider a dissipative 
version of the 0(3) sigma model 

(17) 



ra t M = n J1? [v 2 M + 7f(t,f)] , 



(20) 



where II ^ A 



A — M (M A) for any vector A and the 

magnetization M is subject to the constraint MM = 1. 
The vector noise components are white noises 



< n k (t,x) >= , 

< ?7 fe (*i,^i)V(*2,^2) >= 2fcT L S kl 5(ti - fe) 8(x x - x 2 ) (21) 

for any k, I = 1, 2, 3. 

For a positive topological index n, in the absence of 
noise, the solution of Eq.(|2(]) relaxes to one out of the 
continuous family of degenerate ground states 0] 



M (x,c,a,i,bi) = ( 
(z - ax) ... 



W 



W 



1 + \w\ 2 

(z - a n ) 



w* w 



W* 1 — |V^| 2 



1 + |IT| 2 ' 1 + \W\ 2 



{z-bx)...{z- b„) 



(22) 



parametrized by the complex parameters a's, 6's and c; 
z — x\ — 1X2- There are (An + 2) real parameters. 3 of 
them can be removed by fixing the boundary conditions 
at planar infinity and the remaining global U(l) sym- 
metry [^J. After this reduction there are (An — 1) real 
parameters left, they parametrize a (An — l)-dimensional 
real manifold M n , which is called the moduli space [Q. 
We will denote the real parameters by £ Q with a — 
1, . . . , (An - 1). The static solution M (x, £) has (An - 1) 
independent zero modes dMo(x,£)/dt; a . 

As for kinks, compare Eq.(||), the magnetization can 
be expanded in excited states 



M(t,x) = M [x,Z(t)yi-K 2 + K, 



where 



a#0 



(23) 



(24) 



Note that K ■ Mq = 0. The sigma model analogue of the 
first of Eqs. (0) is 

rwK)^ = 'i a (a) • (25) 

where g a p is a metric tensor on the moduli space 
dM (x^)9M (x,0 



g«m = jdx ^ d ^ 

Time correlations of the projected noises 
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Va(t,0 = / d 2 X ff(t,x) 



can be obtained from the correlations (|2l]) as 

<»?a(*,0 >=0 , 

< V<x(ti,0V(3(t 2 ,0 >= 2kTT g a() (i) 5(ti-t 2 ) 



(27) 



(28) 



The equation ( ]25| ) and the correlations (|2J]) lead to the 
following diffusion equation on the moduli space 



dt 



_d 



_d_ p(t,p 



(29) 



sigma model. Even in the pure sigma model, in which the 
solitons do not interact by potential forces, the thermal 
noise induces entropic forces. In particular, two solitons 
attract each other at large distances and repel at short 
distances, and so the entropic forces tend to bind them 
into bi-solitonic molecules. 
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<?(£) = det[g af3 (Z)}, g^ = {g- 1 )^ and D a = kT/T 
is a diffusion coefficient. P(t,£,) is the probability of 
finding the multisoliton at the point £ on the moduli 
space, such that the average of any function F(£) is 
given by < F >= f Mn d£ P(t,£)F(£). If f are re- 
stricted to a compact area of the moduli space, the prob- 
ability would finally relax to the equilibrium distribu- 
tion P eq (0 ~ y/g{Q- Interpreting this distribution as a 
Boltzman distribution, we can identify the entropic po- 
tential on the moduli space as 



TS($ = -kT In y/g^j 



(30) 



which is not quite unexpected as yj g{£) is a surface ele- 
ment on M n . 

Example. To give a clear but nontrivial example of 
entropic forces, let us restrict our attention to the sub- 
space of Ah defined by W = z ^ A i with real A,B. 
This space represents states of two solitons of the over-g . 2 
all size B. For A > (A < 0) they are located on the 

real (imaginary) axis at the points — A ( — iA). The 
sigma model is scale invariant, so it is not surprising that 
g(A,B) = g(A/B), compare Eqs. (10-12) in Ref. f§, 

0.1 
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VgWBj= ^8[|(A'-|)-(^ y ^tan K )2]sin2 K , 

(31) 

where K = arctan(f?/yl) and K = K[cos(k)], E =q 
E[cos(k)] are complete elliptic integrals of the first and 
second kind respectively. The function ^fg has a unique 
global maximum at A/B w 1.26, compare Fig. 2. For 
a definite scale B, which can be fixed by residual in- 
teractions inducing an effective potential on the moduli 
space Q such as L — S coupling, Zeeman energy and a 
higher derivative (Skyrme) term, the entropic forces tend 
to bind the two solitons Sk into a molecule SA^- At short 
distance there is a logarithmic repulsion. 
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FIGURE 1. The equilibrium diffusion coefficient 
-D(oo) for a kink as a function of temperature kT. Solid 
line - perturbative prediction, see Eq.dlgl); dots - results 
from numerical simulation. 



CONCLUSIONS. 



The numerical simulations confirm our method of ex- 
pansion in powers of VkT, when applied to kinks. We be- 
lieve the same expansion is valid for multisolitons in the 
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FIGURE 2. The probability distribution @ to find 
two skyrmions separated by the distance A as a function 
of A/B. 
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